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MATHEMATICAL  MODELING  of  the  dynamic  response 
of  an  occupant  to  severe  decelerative  forces 
has  become  widespread  in  the  engineering  com- 
munity. This  is  due  to  the  dangers  inherent 
in  human  testing  at  these  levels  as  well  as 
the  inability  to  fabricate  a reliable  dummy 
that  exhibits  motional  fidelity  when  compared 
to  corresponding  human  responses.  Satisfactory 
results  which  would  establish  models  as  a 
reasonable  human  analog  have  not  yet  been 
achieved.  This  limited  success  to  date  is  due 
to  the  difficulties  inherent  in  modeling  the 
complexity  of  the  various  body  interactions  as 
well  as  the  paucity  of  precise  human  data. 
Unacceptable  simulation  results  of  the  head 
and  neck  response  can  usually  be  traced  to  in- 
adequate sophistication  in  modeling  the  head 
and  neck  or  to  errors  propagated  by  deficiencies 
in  modeling  the  thorax  - restraint  system  inter- 
actions. 

Tracking  only  head  motion  during  live 
human  tests  is  insufficient  for  simulation 
validation,  since  one  may  adequately  match 
head  angular  displacement  in  the  inertial  ref- 
erence frame  while  being  inaccurate  in  relative 
head-neck-torso  displacements.  Consequently, 
more  than  just  the  head  segment  should  be  moni- 
tored so  that  relative  positions  of  the  various 
instrumented  segments  can  be  established  and 
the  adequacy  of  the  simulation  ascertained  on 
more  than  single  segment  comparison.  If  one 
assumes  that  the  forces  propagated  to  the 
head  are  transmitted  via  the  first  thoracic 
vertebral  body  (Tl) , then  the  forces  at  T1  and 
not  the  sled  profile  define  the  driving  func- 
tion of  the  head.  Only  when  Tl  is  instrumented 
can  relative  motion  between  the  head  and  torso 
be  determined  with  precision  and  the  simulated 
thoracic-restraint  belt  interaction  verified. 


This  paper  analyzes  data  on  the  dynamic  response 
of  the  living  human  head  and  neck  to  -Gx  (eye- 
balls out)  impact  acceleration  where  the  motion 
of  the  subject's  head  and  Tl  was  monitored  in 
the  midsagittal  plane  with  inertial  instrumen- 
tation and  high  speed  photography  for  confir- 
mation. 

The  Calspan  "3D  Computer  Simulator  of 
Motor  Vehicle  Crash  Victims,"  Ultrasystems 1 
"Crash  Victim  Simulator  - Light  Aircraft"  and 
Boeing  Computer  Services'  "Prometheus"  were 
used  to  provide  estimates  of  occupant  responses. 
These  three  programs  represent  a sample  of  the 
various  complexities  and  approaches  employed 
in  modeling  biodynamic  responses.  Inputs  to 
the  programs  were  made  as  uniform  as  program 
restrictions  allowed  and  outputs  were  compared 
for  consistency  as  well  as  accuracy  in  repli- 
cating the  human  results  obtained. 

METHODS 

INPUT  PARAMETERS  - The  specific  human  runs 
selected  for  possible  simulation  were  all  above 
10  G and  relatively  high  G onsets  (greater  than 
500  G/ second).  It  was  felt  that  at  these  levels 
the  motion  of  the  head  and  neck  would  r.ot  be 
influenced  by  muscular  activity  resulting  from 
anticipation  on  the  part  of  the  test  subjects. 
The  subject  finally  chosen  for  simulation  ex- 
hibited typical  response  characteristics  in 
the  head  angular  and  head  linear  accelerations 
attained  with  the  time  histories  of  the  various 
acceleration  profiles  being  very  similar  to 
those  of  the  group  as  a whole.  Based  on  this 
analysis  the  subject  was  felt  to  be  a represen- 
tative sample  of  the  response  to  be  expected 
at  this  level  of  deceleration. 

Detailed  descriptions  of  the  instrumenta- 
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This  paper  analyzes  data  on  the  dynamic 
response  of  the  living  human  head  and  neck  to 
-Gx  impact  acceleration.  The  Calspan  "3D 
Computer  Simulator  of  a Motor  Vehicle  Crash 
Victim,"  Ultrasystems  "Crash  Victim  Simulator 
Light  Aircraft"  and  Boeing  Computer  Services 
"Prometheus"  were  used  to  provide  estimates 
of  the  responses  monitored.  Inputs  to  the 
programs  were  made  as  comparable  as  program 


restrictions  would  allow.  Outputs  were  com- 
pared to  each  other  as  well  as  to  the  corres- 
ponding human  test  run.  Program  outputs  proved 
to  be  consistent  but  failed  to  adequately 
replicate  human  results.  Inclusion  of  head 
to  neck  articulation  did  not  by  itself  improve 
results.  Relocation  of  the  head  pivot  away 
from  the  occipital  condyles  or  introduction  of 
muscular  activity  was  indicated. 
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tion  and  methods  employed  have  been  previously 
presented- (1,2)*.  Briefly,  both  external 
(whole  body)  and  x-ray  (head  and  neck)  anthro- 
pometric data  were  obtained  for  the  subject 
tested.  The  use  of  x-ray  anthropometry  enabled 
the  fixing  of  the  transducer  mounts  and  asso- 
ciated instrumentation  to  well-defined  anato- 
mical coordinate  systems  (i.e.,  head  and  Tl). 
Locations  of  the  head  C.G.  and  occipital  con- 
dyles were  established  relative  to  the  head 
anatomical  coordinate  system,  permitting  trans- 
formation of  the  transducer  data  to  these  points 
using  rigid  body  mechanics.  Since  both  the 
head  and  spine  were  instrumented,  the  corres- 
ponding head  and  Tl  coordinate  systems  could 
be  tracked  in  time,  enabling  the  head  angular 
displacement  to  be  defined  relative  to  both 
the  inertial  and  Tl  coordinate  systems. 

The  definition  of  the  head  anatomical 
coordinate  system  and  the  head  C.G.  location 
were  based  on  Walker's  analysis,  while  the 
location  of  the  occipital  condyles  was  deter- 
mined from  lateral  anthropometric  x-rays  of 
the  test  subject-(3).  Similarly,  the  mass  and 
mass  moment  of  inertia  of  the  head  and  neck 
were  estimated  from  Walker's  data  and  scaled 
to  the  anthropometric  characteristics  of  the 
subject.  The  neck  was  defined  to  be  a vector 
extending  from  the  Tl  coordinate  system  origin 
to  the  occipital  condyle  point.  Since  both 
the  Tl  origin  and  the  occipital  condyle  point 
could  be  located  in  the  inertial  reference 
system  throughout  the  run,  neck  angle  was  de- 
fined as  the  angle  the  neck  made  with  the  Z 
axis  of  this  system.  Head-neck  (H-N)  angle 
was  taken  to  be  the  angle  the  neck  vector  made 
with  the  head  anatomical  coordinate  system. 
Using  this  approach,  precise  initial  conditions 
could  be  determined  from  photographic  and  x-ray 
data  and  the  validity  of  the  simulation  ascer- 
tained . 

The  mass  and  mass  moment  of  inertia,  link 
lengths  and  C.G.  locations  of  the  other  rele- 
vant body  segments  were  based  on  Dempster's 
work  - (4).  Establishing  joint  locations  in 
the  various  segments  proved  to  be  the  most 
vexing  problem.  The  occipital  condyles  were 
chosen  as  the  articulation  point  between  the 
head  and  neck,  while  the  origin  of  the  Tl 
coordinate  system  was  assumed  to  be  the  pivot 
between  the  neck  and  thorax.  Since  the  occi- 
pital condyles  are  anatomically  distributed 
areas  instead  of  points  and  are  not  in  the  mid- 
sagittal  plane  of  the  head,  the  precise  head 
pivot  location  is  somewhat  inexact  but  was 
taken  to  be  the  occipital  condylar  point  pro- 
jected onto  the  midsagittal  plane  of  the  head. 
The  locations  of  the  other  joints  (sterno- 
clavicular, shoulder,  elbows,  waist,  hips, 
knees,  ankles)  were  estimated  by  scaling 
Dempster's  approximations.  All  joints  were 
assumed  to  be  ball  and  sockets  with  the  excep- 
tion of  the  knees  and  elbows,  which  were  con- 


*  Numbers  in  parentheses  designate  References 
at  the  end  of  paper. 


sidered  to  be  hinged.  Limits  of  joint  motion 
for  the  extremeties  and  the  torso  were  based 
on  data  from  Dempster  as  well  as  Glanville  and 
Kreezer  - (5) . Head  range  of  motion  was  based 
on  Glanville's,  Bowman's  and  Tarriere's  data  - 
(6,  7).  It  should  be  noted  that  Glanville's 
forced  range  of  motion  agrees  favorably  with 
Tarrier's  voluntary  values  whereas  Bowman's 
results  are  more  in  keeping  with  Glanville's 
voluntary  ranges  (Table  1).  The  proportions 
of  the  total  range  accounted  for  by  dorsiflexion 
and  ventrif lexion  are  shown  in  Table  2.  Again, 
there  is  good  agreement  between  Glanville  and 
Tarriere.  The  total  range  of  motion  was  broken 
down  into  neck-torso(N-T)  and  head-neck  (H-N) 
components  using  Bowman's  and  Tarriere's 
analysis  (Table  3).  Ventrif lexion  values  for 
N-T  and  H-N  ranges  are  in  relatively  good  agree- 
ment as  are  dorsiflexion  H-N  values.  Dorsi- 
flexion N-T  ranges,  however,  show  a large 
discrepancy,  with  Bowman's  dorsiflexion  values 
accounting  for  only  20  percent  of  the  total 
N-T  range  of  motion.  The  simulation  values 
used  (Table  4)  were  based  on  Tarriere's  pro- 
portional breakdown  of  the  N-T,  H-N  ranges 
applied  to  the  averages  shown  in  Table  1.  All 
joints  were  considered  to  have  almost  free 
motion  (i.e.,  only  friction  and  viscous  re- 
tarding forces)  up  to  the  joint  limiting  angle, 
where  spring  characteristics  were  introduced 
to  restore  the  segment  excursion  back  within 


Table 

1 - Total  Range 

of  Motion 

Clanvllle 
and  Kreezer 

Bownan  and 
Robbins 

Tarriere 
and  Sapin 

Average 

DORSIFLEXION 

(Backward) 

Forced  77° 

Voluntary  61° 

45° 

ioo® 

83.5° 

VENTRIFLEXION 

(Forward) 

Forced  76° 

Voluntary  60° 

70° 

55° 

80.5° 

Total  Range 

Forced  153° 

Voluntary  121° 

115° 

185° 

169.0° 

Table  2 - Proportions  of  Total  Range  of  Motion 

Glnrrvllle 
and  Kreezer 

Bownan 

Robbins 

Tarriere 

and  Sao In  Average 

DORSIF  LEXION 

Forced  0,5 

Voluntary  0,5 

0.39 

0.54 

0.52 

VERTRIFLEX ION 

Forced  0.5 

Voluntary  0.5 

0.61 

0.46 

0.48 

Table  3 

- Neck  - Torso, 

Neck  - Head  Breakdown 

N-T 

1*  - H 

Bovtr-an  Tarriere 

Robbins  and  Sapln 

Bovtan  Tarriere 

Robbins  and  Sapln 

DORSIFLEXION 

15° 

70° 

30°  30° 

VENTRIFLEXION  60° 

65° 

10°  20’ 

Range 

75° 

135° 

40°  50° 

Table  A - 

Head  and  Neck  Pivot  Joint  Ranges  for  Simulation 

N-T 

DORSIFLEXION  61.95° 

N - H 
26.55° 

Total 

88.5° 

VENTRIFLEXION  61.56° 

18.94° 

80.5° 

Range 

123.51° 

45.49° 

169.0° 
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its  joint  stop  contour.  Estimates  of  the  N-T 
and  H-N  forcing  element  (spring  and  dashpot 
function)  were  obtained  from  Becker's  pre- 
liminary analysis  of  data  from  the  subjects 
tested  - (8) . 

The  dimensions  of  all  body  segment  con- 
tact surfaces  were  based  on  the  subject's 
anthropometry.  Kroell's  data  on  the  impact 
tolerance  of  the  human  thorax  was  used  to 
estimate  the  force  deflection  properties  of 
the  thorax  - (9).  Dynamic  stress-strain 
curves  for  the  restraint  harness  webbing  were 
also  obtained  and  used  as  inputs  to  the  pro- 
gram - (10) . 

3D  Computer  Simulator  of  Motor  Vehicle 
Crash  Victim  (CALS PAM)  - This  program  was 
developed  to  simulate  motor  vehicle  crash 
victims  and  extensive  dummy  validation  efforts 
have  been  conducted  - (11,12).  The  equations  of 
motion  are  formulated  in  a manner  that  allows 
variation  of  the  number  of  segments  and  joints 
without  altering  the  program  structure.  The 
program  is  well  suited  to  study  an  occupant  in 
a high  deceleration  or  crash  environment. 
Interactions  of  the  occupant  with  his  environ- 
ment (i.e.,  seat,  restraint  system,  cabin 
configuration)  are  handled  with  sophistication 
and  this  characteristic  together  with  its  mod- 
ular structure  made  it  a prime  candidate  to 
investigate  the  responses  observed  from  the 
test  subjects. 


The  standard  15  segment,  14  joint  version 
(Figure  1)  was  deemed  inappropriate  in  that 
the  entire  upper  torso  was  represented  as  one 
segment.  Since  a double  shoulder  strap  res- 
traint system  was  used  in  the  runs,  it  would 
seem  probable  that  although  the  shoulders  were 
adequately  restrained,  the  thorax  and  conse- 
quently T1  could  exhibit  greater  excursion  than 
the  shoulders  themselves.  The  configuration 
of  the  upper  torso  was  modified  with  the  inclu- 
sion of  two  sternoclavicular  joints  (ball  arri 
socket)  and  was  broken  down  into  thorax  and 
two  shoulder  elements.  This  modification 
allowed  the  distribution  of  weight  to  be  appor- 
tioned in  a manner  consistent  with  Dempster's 
cadaver  data  and  required  no  source  modifica- 
tion to  the  computer  program.  The  various 
joints  were  located  in  the  relevant  segments 
(using  Dempster's  approximations)  and  this 
16  segment-15  joint  representation  is  shown 
in  Figure  2. 

It  should  be  noted  that  the  program  has 
provisions  for  usage  of  a sophisticated  shoulder 
joint  model  which  basically  consists  of  upper 
torso,  scapula  and  upper  arm  segments  connected 
via  ball  and  socket  joints  and  spring  dampers 
to  simulate  muscles.  Relative  motion  between 
the  scapula  and  upper  torso  segments  is  limited 
by  a fixed  distance  constraint,  simulating  the 
clavicle,  and  a rolling-sliding  constraint 
between  the  respective  contact  surfaces.  Due 
to  the  absence  of  reliable  input  data,  this 
particular  option  of  the  program  was  not 
employed  but  will  be  considered  in  future 
simulations.  Additionally,  a new  numerical 
integrator  (Fleck  Integrator),  which  appears 
to  be  substantially  more  efficient,  has  been 
developed  by  Calspan  but  was  not  used  in  our 
rams  - (13)  . It  is  expected  therefore  that  the 
noted  execution  times  in  this  paper  can  be 
drastically  reduced. 

Crash  Victim  Simulator-Light  Aircraft 
(Ultrasystems)  - The  CVSLA  computer  program 
consists  of  a three-dimensional  model  of  an 
aircraft  seat,  occupant  and  restraint  system 
and  was  developed  as  an  aid  in  the  design  and 
analysis  of  crashworthy  seats  and  restraint 
systems  for  general  aviation  aircraft.  The 
user  has  the  option  of  choosing  among  several 
basic  seat  designs  incorporated  within  the 
program,  ranging  from  a rigid  configuration 
to  one  including  provisions  for  simulating 
plastic  behavior.  A rigid  seat  was  used  for 
this  simulation.  The  modeled  restraint  system 
consisted  of  a lap-belt  and  a double  shoulder 
strap  harness.  Since  this  model  does  not  con- 
sider segment  deformation  (as  is  the  case  for 
the  Calspan),  the  force  elongation  properties 
of  the  shoulder  strap  webbing  were  modified 
to  reflect  the  thoracic  deformation  that  would 
occur  when  these  belts  are  loaded.  The  occu- 
pant is  modeled  via  11  rigid  segments  and  10 
joints  (Figure  3).  In  the  original  version, 
joint  ranges,  segment  dimensions  and  inertial 
properties  were  determined  from  studies  of 
human  body  anthropometry  and  kinematics  - (14) . 
Segment  properties  were  reduced  to  fixed  frac- 
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Fig.  2-16  segment,  15  joint  representation 
used  for  modified  Calspan  occupant  model 


tions  of  total  body  weight  and  stature.  These 
fractions,  along  with  joint  properties  and  a 
multiplication  factor  to  account  for  total 
body  size,  are  incorporated  in  the  program  and 
are  used  to  calculate  the  occupant's  anthropo- 
metric characteristics.  This  approach  greatly 
reduces  the  inputs  to  the  program  although  it 
precludes  the  modeling  of  a specific  test  sub- 
ject's anthropometric  characteristics. 

Consequently,  for  purposes  of  this  paper 
it  became  necessary  to  alter  these  fractional 
calculations  to  insure  that  the  link  lengths, 
mass  distribution,  and  joint  locations  would 
approximate  the  subject  simulated.  This  was 
accomplished  without  affecting  the  operation 
of  the  program,  although  it  circumvented  the 
user  oriented  aspects  of  the  input.  Addition- 
ally, joint  ranges  and  loading  properties  were 
modified  to  agree,  as  close  as  the  various 
methods  of  calculation  would  allow,  to  the 
inputs  used  in  the  Calspan  simulation.  Link 
length  and  C.G.  location  of  the  head  and  neck, 
which  is  modeled  as  a single  rigid  element, 
together  with  the  initial  condition  of  this 
segment  were  matched  as  close  as  possible  to 
those  of  the  human  run.  Outputs  of  the  program 


were  expanded  to  include  the  accelerations  of 
the  neck  link  (Tl)  as  well  as  the  angular 
acceleration  and  velocity  of  the  head  C.G. . 

As  a result  of  these  modifications  the  only 
significant  occupant  differences  from  the 
Calspan  model  was  the  representation  of  the 
head  and  neck  as  a single  segment  (vs.  head 
and  neck  connected  by  a joint)  and  the 
modeling  of  the  upper  torso  as  a rigid  element 
(vs.  thorax  and  two  shoulder  segments  connected 
by  two  sternoclavicular  joints).  This  results 
in  a somewhat  different  distribution  of  masses 
for  the  upper  torso  and  the  head-neck  segment, 
as  well  as  decreased  articulation  between  the 
thorax  and  the  head. 

Prometheus  (Boeing)  - This  program  is  a 
two-dimensional  crash  victim  simulator.  The 
occupant  is  modeled  by  eight  concentrated 
masses  connected  by  seven  links  (Figure  4), 
and  is  restrained  by  lap  and  shoulder  belts. 

The  occupant  interacts  with  the  environment 
through  a non-linear  finite  element  model  of 
the  seat  structure.  The  program  incorporates 
a number  of  user-oriented  features,  and  is 
especially  useful  for  studying  seat  structure 
behavior  during  helicopter  and  fixed-wing 
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a ratio  test  of  0.05  and  a 0.1  ms.  minimum  step 
size,  no  check  on  the  0.1  ratio  test  was  conduc- 
ted. Both  CVSLA  and  Prometheus  calculate  error 
bounds  for  each  variable,  comparing  them  at  each 
time  step  with  the  difference  between  the  pre- 
dicted and  corrected  values.  Calspan  on  the 
other  hand,  has  the  option  of  using  different 
ratio  tests  for  the  various  variables.  If  no 
ratio  test  limit  is  provided  for  a given  variable, 
no  check  is  conducted  on  that  variable  and  the 
convergence  test  for  a given  integration  step  can 
be  passed  despite  possible  inaccuracies  or  even 
instabilities  in  that  quantity.  All  Calspan 
runs  presented  in  this  paper  however,  applied  the 
ratio  test  quoted  to  the  linear  acceleration  of 
the  reference  segment,  as  well  as  to  the  angular 
acceleration  components  of  all  16  segments  mod- 
eled. Intra-program  outputs  for  the  various 
ratio  tests  were  compared  for  consistency  and 
Figures  5,  6,  and  7 show  the  results  for  the  head 
angular  acceleration,  one  of  the  most  sensitive 
output  variables.  All  Calspan  and  Prometheus 
outputs  were  virtually  unchanged  over  the  differ- 
ent ratio  limits  while  CVSLA  showed  slight  varia- 
tion for  the  0.05  ratio  test,  but  appeared  sta- 
ble. The  minimum  permissible  step  sizes  for 
Calspan  and  Prometheus  were  never  reached,  indi- 
cating that  the  variables  passed  the  ratio  tests 
without  having  to  resort  to  the  smallest  step 
size  provided.  The  CVSLA  output  was  not  exten- 
sive enough  to  determine  whether  the  minimum  al- 
lowable time  intervals  were  approached,  but  the 
appreciable  increases  in  program  execution  times 
with  the  lower  ratio  tests  did  preclude  the  chro- 
nic use  of  the  minimum  step  size.  Additionally, 
the  virtually  identical  outputs  for  the  0.1  and 
0.01  ratio  tests  (using  minimum  step  sizes  of  0.1 
and  0.01  ms.  respectively)  indicates  that  the 
minimum  time  interval  for  these  two  runs  was  ade- 
quate. Inter-program  comparisons  presented  below 
used  ratio  test  limits  of  0.05,  0.1,  and  0.1  for 
Calspan,  CVSLA,  and  Prometheus,  respectively. 

The  CPU  time  necessary  to  execute  the  programs 
depended  strongly  on  the  integration  parameters. 
With  all  three  programs  integrating  to  250  ms. 
typical  CPU  execution  times  on  a CDC  6600,  for 
the  ratio  test  limits  shown,  ranged  from  296  to 
754  sec.  for  Calspan,  from  132  to  354  sec.  for 
CVSLA,  and  from  36  to  49  sec.  for  Prometheus. 

Comparison  Criteria  - With  the  exception  of 
minor  program  modifications  to  CVSLA  to  make 
inputs  comparable,  no  rewriting  of  source  code 
wa3  conducted  in  any  of  the  programs  to  change 
the  respective  mechanical  idealizations.  The 
programs,  as  delivered,  were  thoroughly  checked 
out;  sample  problems  were  run,  and  the  results 
compared  to  those  obtained  on  other  computer 
systems.  Since  results  proved  to  be  consistent, 
the  programs  were  assumed  to  be  error  free  and 
the  differences  in  program  results  were  attri- 
buted to  the  respective  mechanical  idealizations 
employed.  Although  convergence  for  the  varia- 
bles was  checked,  possible  program  errors,  re- 
sulting in  inaccuracies  of  the  numerical  solu- 
tions, cannot  be  absolutely  ruled  out  as  contri- 
butors to  the  differences  in  results.  To  elimi- 


nate all  possibility  of  program  errors,  all 
three  simulators  would  have  to  be  run  and  com- 
pared under  identical  conditions  and  mechanical 
idealizations  (same  number  cf  segments,  identi- 
cal loading  and  unloading  functions,  etc.).  The 
only  major  differences  then  would  be  the  various 
integration  procedures  employed.  A direct  com- 
parison of  results  under  these  conditions  would 
isolate  possible  programming  deficiencies  a3 
well  as  make  a comparison  of  execution  times 
more  meaningful,  since  identical  idealizations 
and  comparable  accuracy  would  have  been  used. 

Due  to  the  user  oriented  nature  of  this  paper, 
this  analytic  approach  was  not  taken. 

RESULTS 

In  order  to  make  a direct  comparison  of  the 
three  programs  possible,  the  head  pivot  and  cla- 
vicular joints  in  the  Calspan  model  were  locked. 
Since  Prometheus  does  not  allow  for  a footrest, 
the  modeled  environment  consisted  solely  of  seat 
back,  seat  pan  and  floor  surfaces.  Consequently, 
with  the  exception  of  a somewhat  different  dis- 
tribution of  masses,  the  three  representations 
of  the  test  subject  and  interacting  surfaces  were 
as  uniform  as  the  various  programs  would  allow. 
Differences  in  results  should  be  due  primarily  to 
the  methodologies  used  in  modeling  occupant-res- 
traint interactions  and  joint  loading  and  un- 
loading characteristics.  The  results  proved 
to  be  relatively  consistent,  although  signifi- 
cant differences  were  evident,  especially  in 
the  predicted  belt  loads.  For  the  lap  belt 
load  (Figure  8)  both  CVSLA  and  Prometheus  dis- 
played higher  amplitudes  than  Calspan,  although 
the  loading  characteristics  were  quite  similar. 
Calspan  and  Prometheus  reached  their  peak  forces 
at  approximately  the  same  time  (80  ms.),  where- 
as the  corresponding  CVSLA  value  was  delayed 
by  35  ms.  This  delay  was  probably  due  to  motion 
of  the  occupant  under  his  restraint,  since  both 
Calspan  and  Prometheus  assumed  the  lap  belt  to 
be  anchored  to  the  lower  torso,  while  CVSLA 
allowed  for  submarining.  Unloading  character- 
istics of  the  lap  belt  were  similar  for  CVSLA 
and  Prometheus  whereas  Calspan  had  a much  more 
gradual  decrease  in  forces  generated. 

Shoulder  strap  loads  (Figure  9)  also  demon- 
strated differences  in  occupant  restraint  inter- 
actions. CVSLA  and  Prometheus  had  similar  fre- 
quencies of  oscillation  with  the  latter  damping 
out  much  more  rapidly.  This  damping  is  also 
evident  in  Calspan,  although  the  frequency  was 
higher  than  in  the  other  two,  indicating  a 
stiffer  interface.  It  must  be  remembered  that 
these  variations  are  primarily  due  to  the  dif- 
ferent force  deflection  models  employed  by  the 
respective  programs.  Although  all  three  pro- 
grams use  tabular  input  in  defining  loading 
characteristics,  they  differ  in  their  treatment 
of  the  unloading  algorithm.  Prometheus  employs 
a fixed  tabular  input  and  CVSLA  a fixed  poly- 
nomial approximation.  Calspan  on  the  other 
hand  uses  a dynamic  unloading  process  based  on 
energy  absorption  and  energy  restitution  coef- 
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ficients.  It  appears  that  the  bulk  of  the  vari- 
ation in  results  can  be  attributed  to  the  dif- 
ferent treatments  of  hysteresis. 

This  variation  in  loaning  of  the  shoulder 
straps  was  reflected  in  the  T1  accelerations 
(Figure  10).  CVSLA  and  Calspan  agreed  well  in 
initial  accelerations  achieved,  and  time  of 
peak  loading,  but  different  frequencies  of  os- 
cillation were  again  evident,  with  the  latter 
being  more  damped.  Prometheus  proved  to  be 
similar  to  Calspan  but  indicated  a greater 
initial  G level  than  was  predicted  by  the  other 
two. 

Since  the  head-neck  system  was  rigid  in 
all  cases,  the  neck  angular  displacement  can 
be  considered  coincident  with  head  angular  dis- 
placement. Neck  angular  displacement  results 
(in  the  inertial  reference  frame)  are  shown  in 
Figure  11.  Since  neck-thorax  joint  ranges  were 
identical,  variability  in  neck  angular  displace- 
ments attained  were  a result  of  differences  in 
the  respective  forcing  functions  (T1  accelera- 
tion) as  well  as  variability  in  the  joint  load- 
ing calculations  once  the  limits  had  been 
reached.  Calspan  and  Prometheus  results  dis- 
played a larger  excursion  than  CVSLA,  indica- 
ting a softer  joint  stop. 

Head  angular  velocity  and  acceleration 
(Figures  12  and  13)  were  similar  in  both  the 
magnitudes  attained  as  well  as  in  their  time 
histories.  Calspan  and  CVSLA  head  angular  velo- 
cities were  comparable  in  the  initial  stages 
with  differences  in  the  later  values  probably 
due  to  unloading  of  the  joint  stops.  Prometheus 
although  attaining  a larger  initial  value,  gen- 
erally replicated  the  results.  Trends  in  head 
angular  accelerations  were  also  consistent  in 
all  the  simulations.  All  exhibited  the  charac- 
teristic initial  bi-phasic  time  history  followed 
by  a reversal  in  head-neck  rotation  as  the  neck 
joint  unloaded.  Oscillations  in  linear  accel- 
eration evident  at  T1  were  relatively  damped 
out  upon  reaching  the  head  C.G.  (Figure  14) . 
Prometheus  again  tended  to  give  the  highest  pre- 
dicted G level  attained  but  the  trends  exhibited 
by  all  three  programs  were  similar. 

Although  results  described  do  indicate 
good  correlation  in  trends  and  magnitudes,  the 
relevant  question  as  to  the  validity  of  these 
simulations  still  had  to  be  addressed.  Pre- 
dicted neck  angular  displacements  were  compared 
to  the  human  data  (NAMRL)  and  of  the  three  pro- 
grams Calspan  results  (L0CK-HP,CL)  were  found 
to  give  the  closest  approximation  (Figure  15). 

It  is  crucial  to  remember  that  neck  angle  was 
defined  as  the  angle  a vector  extending  from  T1 
origin  to  the  occipital  condylar  point  makes 
with  the  inertial  reference  frame.  What  this 
result  indicates  is  that  the  simulation  was 
quite  accurate  in  predicting  the  relative  dis- 
placement between  T1  and  the  occipital  condylar 
point.  Consequently,  if  the  head-neck  system 
were  in  fact  rigid,  the  simulation  (as  far  as 
predicting  head  and  neck  motion)  could  be  con- 
sidered a success  in  that  the  head  angular  dis- 
placement would  by  definition  be  just  as  accu- 


rate. Human  head  to  neck  angular  displacement 
(that  is,  the  angle  the  neck  vector  makes  with 
the  head  anatomical  coordinate  system)  is  shown 
in  Figure  16.  This  analysis  showed  that  by  the 
definition  of  head  and  neck  used  in  this  paper 
there  did  exist  head-neck  dorsif lexion,  reaching 
a value  of  approximately  37  degrees  (NAMRL) . 
Consequently,  since  relative  head-neck  motion 
is  evident,  articulation  between  the  head  and 
neck  segments  must  be  provided.  Since  only 
Calspan  had  this  facility,  both  the  clavicles 
and  the  head  pivot  were  sequentially  unlocked 
and  the  program  was  rerun  for  these  additional 
degrees  of  freedom  (LOCK-HP  and  UNLOCKED 
respectively).  The  neck  angular  displacement 
results  (Figure  15)  were  not  significantly 
changed  from  the  locked  case.  The  absolute 
head  angular  displacement  (Figure  17)  preceded 
and  exceeded  that  of  the  test  case  and  since 
the  neck  angles  agreed  well,  this  error  lies 
in  the  relative  head  to  neck  motion.  The 
simulation  head-neck  angle  did  not  achieve 
the  degree  of  dorsiflexion  that  was  indicated 
in  the  test  data  and  demonstrated  appreciable 
head-neck  ventriflexion  once  the  neck  torso 
joint  stop  began  loading  (Figure  16). 

It  should  be  noted  that  only  minimal 
head  to  neck  dorsiflexion  was  attained  by  un- 
locking of  the  head  pivot.  In  fact,  for  the 
first  140  ms.  the  absence  of  significant  head- 
neck  motion  would  indicate  that  UNLOCKED  re- 
sults should  be  in  good  agreement  with  the 
two  locked  joint  simulations.  Head  angular 
velocity  and  acceleration  did  not  differ  sig- 
nificantly for  the  locked  vs.  unlocked  cases 
(Figures  18  and  19),  although  unlocking  of 
the  head  pivot  and  clavicles  did  delay  the 
corresponding  curve.  In  relation  to  the  test 
data,  the  head  angular  velocity  of  the  unlocked 
joints  case  was  of  lower  magnitude  with  slower 
rise  and  decay  times  while  the  head  angular 
accelerations  failed  to  replicate  the  high 
initial  amplitude  and  following  overshoot. 
Results  of  the  head  resultant  linear  acceler- 
ation were  generally  poor  in  that  although  the 
bi-phasic  aspect  of  the  test  curve  was  evident, 
the  magnitudes  differed  significantly  (Figure 

20) .  Tl  linear  acceleration  proved  reasonable 
for  the  first  half  of  the  simulation,  exhibiting 
similar  trends  and  amplitudes  achieved  (Figure 

21) . 

SUMMARY  AND  CONCLUSIONS 

Although  simulation  results  from  the 
various  programs  agreed  well,  their  accuracy 
in  replicating  the  test  data  was  generally 
poor.  Inter -program  correlations  could  be 
improved  by  further  adjustments  to  the  input 
parameters  but  the  basic  problem  of  inconsis- 
tency with  test  results  would  still  remain  un- 
resolved. Possible  improvements  to  increase 
sensitivity  include  greater  sophistication  in 
modeling  occupant-environment  (seating  and 
restraint  system)  interactions.  Submarining 
and  contacts  with  cabin  surfaces  affect 
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results  and  should  consequently  be  included 
in  models  that  hope  to  address  general  appli- 
cations. Although  user  oriented  features  have 
some  very  obvious  advantages,  the  loss  of  the 
ability  to  model  specifics  can  decrement  usage 
as  well  as  degrade  the  results  obtained.  The 
intended  usage  of  the  programs  will  affect  the 
degree  of  sophistication  required.  If  struc- 
tural response  is  of  prime  interest,  then  occu- 
pant simulation  need  not  approach  the  comp- 
lexity that  would  otherwise  be  warranted.  If 
that  is  the  case  however,  conclusions  about 
the  response  of  that  occupant  to  decelerative 
pulses  should  be  made  with  great  caution. 

Lack  of  precise  input  parameters  still 
represents  one  of  the  most  vexing  problems. 
Relatively  small  errors  in  joint  ranges  and 
loading  and  unloading  properties  as  well  as 
force  deformation  characteristics  of  the 
various  interacting  surfaces  can  significantly 
change  results.  Since  great  care  was  taken 
in  making  inputs  to  the  three  programs  com- 
parable, results  were  generally  consistent 
if  not  accurate.  Additional  sophistication 
is  warranted  since  the  resulting  simulation 
results  did  not  approximate  those  of  the  test. 
The  inclusion  of  additional  articulation  be- 
tween the  head  and  neck  did  not  resolve  the 
problem  in  that  greater  accuracy  was  not 
achieved.  This  is  not  to  say  that  simulation 
of  relative  motion  between  the  head  and  neck 
segments  is  unnecessary.  Human  head-neck 
values  clearly  indicated  that  relative  motion 
between  the  head  and  neck  does  in  fact  take 
place.  Since  increased  articulation,  by  it- 
self, did  not  provide  the  expected  displace- 
ments, other  modifications  have  to  be  con- 
sidered. It  must  be  remembered  that  neck 
angular  displacements  agreed  very  well  and 
consequently  discrepancies  in  the  head  angular 
displacement  must  be  accounted  for  relative  to 
the  occipital  condyles.  Although  the  occipital 


condyles  are  generally  considered  the  anatomical 
pivot  point,  it  need  not  necessarily  be  the 
mathematical  one.  Since  the  head  does  not 
rotate  on  the  neck  at  a well  defined  point  but 
slides  and  rolls  on  the  atlas,  the  center  of 
rotation  of  this  complex  motion  would  define 
the  precise  pivot.  If  this  theoretical  pivot 
lies  somewhere  above  the  head  C.G.,  then  the 
demonstrated  head-neck  dorsiflexion  could  be 
achieved  due  to  the  revised  geometry.  Addi- 
tionally, the  representation  of  the  neck  as  a 
rigid  segment  is  an  oversimplification  that 
would  affect  head  results.  The  introduction 
of  a flexible  neck  would  probably  elevate 
(due  to  whiplash)  the  acceleration  levels 
achieved  but  should  not  in  itself  reverse 
direction  of  motion. 

The  time  history  of  the  human  data  showed 
that  peak  Tl  linear  acceleration  was  followed 
by  peak  angular  and  linear  accelerations  of 
the  head,  indicating  that  forces  to  the  head 
are  propagated  through  Tl.  It  is  interesting 
to  note  that  maximum  head-neck  dorsiflexion 
is  attained  at  approximately  85  ms.,  and  that 
peak  head  linear  and  angular  acceleration 
occur  at  75  and  79  ms.,  respectively.  Since 
the  peak  accelerations  and  maximum  head-neck 
angular  displacement  (dorsiflexion)  are  almost 
coincident  in  time,  these  peak  accelerations 
appear  to  be  caused  by  head-neck  pivot  loading. 
This  theorized  loading-unloading  of  the  head- 
neck  joint  is  consistent  with  head  angular 
acceleration  data  where  unloading  produces 
the  overshoot  exhibited  at  approximately  110 
ms.  Additionally,  a sharp  spike  is  also  evi- 
dent in  the  Tl  linear  acceleration  profile 
at  approximately  80  ms.,  further  substantiating 
possible  head-neck  loading.  Since  variability 
between  simulation  (UNLOCKED)  and  human  test 
results  became  evident  approximately  40-60  ms. 
into  the  pulse,  muscular  contraction  (whose 
response  time  would  be  in  this  domain)  could 
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account  for  the  dorsiflexton  of  the  head  rela- 
tive to  the  neck.  Alterations  to  the  program 
to  substantiate  this  hypothesis,  together  with 
inclusion  of  a flexible  neck,  will  be  under- 
taken before  sensitivity  analysis  and  parameter 
adjustments  will  be  conducted. 
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